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We report ab-initio calculations for the electronic structure of organic charge transfer salts k- 
(ET) 2 Cu[N(CN) 2 ]Br, k-(ET) 2 Cu[N(CN) 2 ]I, k"-(ET)2Cu[N(CN) 2]C1 and k-(ET) 2 Cu 2 (CN) 3 . These 
materials show an ordering of the relative orientation of terminal ethylene groups in the BEDT- 
TTF molecules at finite temperature and our calculations correctly predict the experimentally 
observed ground state molecular conformations (eclipsed or staggered). Further, it was recently 
demonstrated that the ethylene endgroup relative orientations can be used to reversibly tune k- 
(ET) 2 Cu[N(CN) 2 ]Br through a metal-insulator transition. Using a tight-binding analysis, we show 
that the molecular conformations of ethylene endgroups are intimately connected to the electronic 
structure and significantly influence hopping and Hubbard repulsion parameters. Our results place k- 
(ET) 2 Cu[N(CN) 2 ]Br in eclipsed and staggered configurations on opposite sides of the metal-insulator 
transition. 

PACS numbers: 71.10.Fd, 71.15.Mb, 71.20.Rv, 74.70.Kn 


Quasi-two dimensional charge transfer salts k- 
(BEDT-TTF) 2 V, where BEDT-TTF stands for bis- 
ethylenedithio-tetrathiafulvalene, often abbreviated as 
ET, constitute a fascinating family of materials due to 
their rich phase diagrams comprising metallic, supercon¬ 
ducting, Mott insulating and spin-liquid phases mi- 

These electronic properties of ET-based materials 
are very sensitive to disorder. Irradiation experiments 
have shown that lattice disorder lowers the of k- 
(ET) 2 Cu(SCN )2 [S] and causes electron localization in 
K-(ET) 2 Cu[N(CN) 2 ]Br [S]. Early on it was also real¬ 
ized that ET molecules have intramolecular degrees of 
freedom, namely the configurations of their two ethylene 
endgroups (see Fig. |^, which can either be aligned par¬ 
allel (eclipsed, E) or canted (staggered, S) [7HT7]. The 
energetically favorable configuration is not universal for 
different anions X and packing motifs. For some mate¬ 
rials a glass-like freezing of the ethylene endgroups upon 
cooling has been observed [THHHEHini- 

Especially the first ambient pressure ET-based super¬ 
conductor /3-(ET)2l3 attracted a lot of interest, because 
its Tc can be enhanced from 1.5 K to 8 K by forcing the 
ET molecules, which are endgroup disordered at ambient 
pressure, to assume staggered configuration through ap¬ 
plication of shear and pressure [7Hni uni mi- Recently, 
it was shown that ethylene endgroup disorder can be 
used to reversibly tune K-(ET) 2 Cu[N(CN) 2 ]Br through 
a metal to insulator transition [ITj [22l [23] . 

It is believed that materials k-(ET) 2 X have a common 
phase diagram, which is mainly controlled by the value 
of the on-site Coulomb repulsion U over the electronic 
bandwidth |24] . Changes in physical properties in the 
presence of ethylene endgroup disorder have so far been 
interpreted as a consequence of lattice disorder, with the 
exception of recent scanning tunneling spectroscopy ex¬ 
periments [25] . Surprisingly, the effect of different ethy¬ 


lene endgroup configurations on the electronic structure, 
and especially the electronic bandwidth, of k-(ET) 2 V 
has only been investigated in a single material using the 
extended Hiickel method m, while calculations for ET 
molecules and dimers in vacuum are available HZIIISIIIT]. 
A preceding study for «;-(ET) 2 Cu[N(CN) 2 ]Br focused on 
energetics rather than bandstructure effects [l6] . 

Here we examine the electronic structure of endgroup 
ordered crystals in both possible configurations for vari¬ 
ous members of the k-(ET) 2 X family of materials using 
ah-initio density functional theory (DFT) calculations. 
We construct a 3/4-filled low-energy effective Hamilto¬ 
nian using a Wannier downfolding scheme and relate the 
resulting model parameters to the relative orientation of 
terminal ethylene groups. Our results show that ethylene 
endgroup configurations of ET molecules influence the 
electronic bandwidth of all resulting crystalline materials 
investigated here. Finally, we offer a simplified interpre¬ 
tation of our results in terms of the Hubbard model on the 
anisotropic triangular lattice and relate our findings to re¬ 
cent experiments, especially pointing out the possibility 
that changes in the electronic structure through ethylene 
endgroup disorder and strongly enhanced electron corre¬ 
lation are relevant in addition to commonly considered 
lattice disorder. 

We used ab-initio density functional theory (DFT) 
calculations within an all-electron full-potential local 
orbital (FPLO) |28| basis to calculate the electronic 
bandstructure. For the exchange-correlation functional 
we employed the generalized gradient approximation 
(CCA) [^. All calculations were converged on 6 x 6 x 6 
fc-point grids. 

Tight-binding models were obtained from projective 
molecular orbital Wannier functions as implemented in 
FPLO [30]. We have shown previously that, for the 
materials of interest here, this method yields near per- 
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feet representations of the low-energy DFT bandstruc- 
ture and avoids ambiguities from band fitting proce¬ 
dures [21] . The resulting tight-binding Hamiltonian H = 

3/4-filled and consists of four 
bands corresponding to the four ET sites indexed by i,j. 
Each site is located at the center of the inner C-C bond 
of an ET molecule. Subsequently, we relate those ab- 
initio calculated parameters with an effective two-band 
half-filled Hubbard model on the anisotropic triangular 
lattice. 

Experimental crystal structures including both possi¬ 
ble ethylene endgroup configurations are available m 
[31[35j for k-(ET) 2 Cu[N(CN) 2 ]I, k"-(ET)2Cu[N(CN) 2]C1 
and k-(ET) 2 Cu 2 (CN )3 (T = 200 K), while for k- 
(ET) 2 Cu[N(CN) 2 ]Br only one endgroup configuration 
has been reported in the literature [T2|. For the latter 
case, we set up the missing ethylene endgroup orientation 
by hand. For k-(ET) 2 Cu[N(CN) 2 ]I we used the crystal 
structure determined at T = 127 K if not denoted other¬ 
wise. The structure at T = 295 K is used as a consistency 
check. 

In Figure we show as an example the displacement 
ellipsoids resulting from the structure determination of 
k"-(ET)2Cu[N(CN) 2]C1. These represent the uncertainty 
of the experimentally determined structure with respect 
to the position of every individual atom. Displacements 
on the endgroup carbon atoms are largest, but the un¬ 
certainty of the sulphur atom positions next to the ethy¬ 
lene endgroups is also significant. Therefore, we chose 
to relax the entire ethylene endgroups together with the 
two neighboring sulphur atoms in both staggered and 
eclipsed endgroup configurations for all materials inves¬ 
tigated here. For comparison of the relaxed crystal struc¬ 
tures with eclipsed and staggered endgroups see Fig. 
Other atomic positions and the unit cell parameters were 
left untouched. Note that hydrogen positions belonging 
to the ethylene groups were not measured in experiment 
and had to be inserted manually. 

For the structural relaxation we used the projector 
augmented wave (PAW) method [34] as implemented in 
GPAW [3S]. We optimized the terminal ethylene end- 
groups using 2x2x2 fc-point grids and GGA exchange- 
correlation functional until forces were below 10 meV/A. 

It turns out that both endgroup configurations locally 
minimize the total DFT energy, so that stable crys¬ 
tal structures can be obtained for both staggered and 
eclipsed orientations from a DFT relaxation. 

It is experimentally known that different charge trans¬ 
fer salts can either have staggered or eclipsed endgroups 
as their lowest energy configuration. To confirm the va¬ 
lidity of our relaxed structures, we first calculate the en¬ 
ergy difference /S.E between staggered and eclipsed us¬ 
ing FPLO. The results are displayed in Table [l| Val¬ 
ues for the energy differences calculated with GPAW 
are in good agreement. The energy ordering of stag¬ 
gered and eclipsed configurations comes out correctly 
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FIG. 1. (Color online) ET molecule with displacement ellip¬ 
soids from the experimental structnre determination of k ”- 
(ET) 2 Cn[N(CN) 2 ]Cl. The size of the ellipsoids scales with 
the amount of uncertainty regarding the respective atomic 
position. Ethylene endgroup disorder is included on the right 
end of the molecule and given by fractionally occupied posi¬ 
tions. The bars above the molecule indicate which parts of 
the ET molecules are relaxed in our DFT calculations and 
which atomic positions remain at their experimental values. 


TABLE I. Lowest energy conhgurations of the ethylene end- 
groups determined from our DFT relaxed structures. The 
energy difference /S.E to the high energy configuration is cal¬ 
culated per ET molecule from the EPLO total energies. 



configuration 

AA in meV 

K-(ET)2CU2(CN)3 

staggered 

130 

k"-(ET)2Cu[N(CN)2]C1 

eclipsed 

72 

Ac-(ET) 2 Cu[N(CN) 2 ]Br 

eclipsed 

110 

«:-(ET) 2 Cu[N(CN) 2 lI 

staggered 

38 


for all investigated materials with the exception of k- 
(ET) 2 Cu[N(CN) 2 ]I at 295K, where the energy differ¬ 
ence is the smallest and the distribution of endgroup 
configurations measured experimentally was found to be 
51:49 [m [311133|. The value for AE we determined is 
the energy difference between the two local minima of the 
energy corresponding to staggered and eclipsed configu¬ 
rations, which is not to be confused with the activation 
energy. The latter denotes the height of the potential 
barrier between those minima, which can be significantly 
larger than AE m US- Instead, our DFT calculated 
values constitute a lower bound for the activation en¬ 
ergy: k-(ET) 2 Cu[N(CN) 2 ]I with the smallest energy dif¬ 
ference is known to be completely endgroup disordered at 
room temperature [H] and k"-(ET) 2 Cu[N(CN) 2 ]C 1 con¬ 
tains about 20% disorder at room temperature [31] , while 
in the other two materials the amount of endgroup dis¬ 
order is in a range of few percent. 

The electronic bandstructures obtained from the 
molecular Wannier function analysis of the DFT results 
are shown in Fig. for both eclipsed and staggered 
ethylene endgroup configurations. In all bandstructures 
shown, the difference between staggered and eclipsed con¬ 
figurations lies in the electronic bandwidth. Going from 
eclipsed to staggered, the overall bandwidth increases, 
while the width of the two bands closest to the Fermi 
level decreases. 
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FIG. 2 . (Color online) Electronic bandstructure of all investigated materials for staggered and eclipsed ethylene endgroup 
configurations. Staggered ethylene endgroups produce a larger overall bandwidth than eclipsed ones, but reduce the width of 
the two bands closest to the Fermi level. 


We analyse these bandstructures using a minimal 
model [36j for a K-packed layer of individual ET 
molecules (Fig. using the four largest parameters 
also commonly denoted as {bi,p,b 2 ,q). Note 
that our tight binding Hamiltonians also include small 
longer range hoppings. The four largest hopping param¬ 
eters are given in Tab.jn] Parameters ti ( 6 i) and ts ( 62 ) 
decrease from staggered to eclipsed configurations, while 
t 4 (q) increases and ^2 (p) remains about constant. 

The molecular Wannier functions for an ET molecule 
in eclipsed and staggered configuration are shown in 
Fig.i Although the molecular Wannier function hardly 
resides on the terminal ethylene groups, overlaps with 
neighboring ET molecules are influenced by the config¬ 
uration of the endgroups through the direction of their 
bonds with the neighboring sulphur atoms. 

Especially the hopping is strongly enhanced, be¬ 
cause the tails of the Wannier function on the neigh¬ 
boring ET molecule are enlarged. The Wannier func¬ 
tions remain largely unaltered in the direction that cor¬ 
responds to t 2 . Therefore, this parameter largely remains 
constant. In the direction of t 4 the tail on the neighbor¬ 
ing ET molecule is enhanced in the staggered configura¬ 
tion, but the Wannier function on the central molecule 
turns away from this tail because of the altered suplhur- 
ethylene bond direction. Consequently, t 4 is reduced. 
The relative changes in ti are rather small, which is con¬ 
sistent with our analysis of the Wannier functions. 

The effect of the altered hopping integrals and 
its connection to the metal-insulator transition in k- 
(ET) 2 Cu[N(CN) 2 ]Br can also be understood in terms of 



FIG. 3 . (Color online) Molecular Wannier function of an 
ET molecule in staggered and eclipsed configuration for k- 
(ET)2Cu[N(CN)2]Br. The arrows denote the directions of 
dominant hopping processes in the individual molecule model 
with {ti ,t2 jts,14). 
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TABLE II. Values of the molecule model parameters {t\, 12 , 13 , 14 ) in meV, also commonly denoted as ( 6 i,p,fe 2 ,g). Dimer model 
parameters are given as ratios t' jt and U/t calculated from (ti,f 2 ,t 3 ,t 4 ) using formulas described in the text. The second column 
states the experimental ground state of the respective material (I=insulator, M=metal, SC=superconductor). The x in the 
fourth column marks the low energy configuration of the ethylene endgroups. 






tl 

t 2 

ts 

t4 

t'/t 

U/t 

K-(ET)2Cu2(CN)3 

I 

eclipsed 


167 

84.9 

70.4 

30.3 

0.61 

5.8 



staggered 

X 

176 

78.0 

81.4 

18.7 

0.84 

7.3 

«;"-(ET)2Cu[N(CN)2]C1 

M 

eclipsed 

X 

174 

97.3 

50.5 

35.9 

0.38 

5.2 



staggered 


188 

93.4 

64.0 

26.6 

0.53 

6.3 

«-(ET) 2 Cu[N(CN) 2 ]Br 

SC 

eclipsed 

X 

178 

99.0 

59.5 

35.8 

0.44 

5.3 



staggered 


187 

97.1 

70.2 

24.9 

0.58 

6.1 

k-(ET) 2 Cu[N(CN) 2 ]I, 127 K 

M 

eclipsed 


152 

101 

47.2 

29.2 

0.36 

4.7 



staggered 

X 

170 

99.0 

52.4 

18.9 

0.44 

5.8 

«:-(ET) 2 Cu[N(CN) 2 ]I, 295 K 

M 

eclipsed 


153 

92.0 

49.9 

31.5 

0.40 

4.9 



staggered 


164 

92.0 

54.8 

22.2 

0.48 

5.8 



FIG. 4. (Color online) Parameter ratios t'/t and U/t in the 
effective dimer model. The direction of the arrows indicates 
the direction of change in the model parameters going from 
the low to the high energy configuration of the ethylene end- 
groups. Eclipsed (E) endgroup configurations are located on 
the left side (small U/t) of the plot, while staggered (S) end- 
group configurations are located on the right side (large U/t). 

the simplified two-band dimer model. We calculate the 
parameters t, t' and U of the corresponding anisotropic 
triangular lattice Hubbard from the usual geometric for¬ 
mulas t = (|t 2 | + \t 4 \)/ 2 , t' = \H\/2 and use the rough 
estimate U « 2|ti| for the intradimer Hubbard repul¬ 
sion [HHl HTT] . 

Recently, there has been some criticism of this method 
of obtaining a t, t' and U dimer model, because the esti¬ 
mate for U relies on the fact that the unscreened inter- 
molecular Coulomb repulsion Vi within an (ET )2 dimer 
vanishes. In DFT parametrizations of this model Vi was 
however found to be non-negligible |26l EH EH]. 

Calculations based on the constrained random phase 


approximation additionally showed sizeable screening ef¬ 
fects [39] in the dimer approximated model not accounted 
for in our method. In iron-based materials it is however 
a well-known problem that different, even very sophisti¬ 
cated, approximations made in the determination process 
can lead to wildly different parameter sets [lOHlH]. 

Recent Wannier function analysis anyway points to 
the fact that a 3/4-filled Hubbard model of individual 
ET-molecules is more appropriate to understand fine de¬ 
tails of k-ET materials Eiiiia. Unfortunately, there are 
only few examples of many-body calculations based on 
such multi-band Hamiltonias |44l [45]. Alternatively, a 
half-filled extended Hubbard model including non-local 
Coulomb repulsion and exchange interactions has been 
proposed |i5] . 

We however emphasize that so far dimer approxi¬ 
mated Hamiltonians and subsequent many-body calcu¬ 
lations showed remarkable success in explaining some of 
the qualitative properties of ET based materials [13 [48] 
and should be sufficient for the case we investigate here. 
Note that our estimates agree especially well with more 
elaborate calculations including screening effects EU [33] ■ 
Newly discovered effects like multiferroic behavior [50] 
or the still unsettled nature of superconductivity [45l l5T] 
however certainly call for more realistic approaches. 

Fig. a shows the result of our dimer model estimate. 
The change from eclipsed to staggered ethylene group 
configuration universally increases both the frustration 
t'/t and the relative strength of the Hubbard repulsion 
U/t, i.e. U over the bandwidth. 

Comparison of our findings with cellular dynamical 
mean-field theory m and exact diagonalization results 
for the anisotropic triangular lattice Hubbard model |3H] 
explains why K-(ET) 2 Cu[N(CN) 2 ]Br can be tuned into 
a Mott insulating state by activating [22l [23] the en¬ 
ergetically less favorable endgroup configuration: First, 
the material in its lowest energy configuration is already 
close to a Mott insulating phase. Second, the lowest 
energy configuration is the eclipsed one, so U/t can be 













5 


strongly increased by activating the staggered configura¬ 
tion. Therefore, the system crosses the phase transition 
line and a Mott insulator is realized. 

Note that the phase diagram of the anisotropic trian¬ 
gular lattice Hubbard model is not entirely settled and 
slightly different results have been obtained using other 
numerical methods [5H [53] . 

In conclusion, we demonstrated that DFT reliably re¬ 
produces the ground state ethylene endgroup configura¬ 
tion for various K-phase materials. While previous discus¬ 
sion of endgroup conformations in the literature consid¬ 
ered only lattice disorder a relevant issue, we have shown 
that the relative orientation of ethylene endgroups within 
ET molecules crucially influences the electronic band¬ 
width of K-type organic charge transfer salts. Switch¬ 
ing an ET molecule from eclipsed to staggered config¬ 
uration decreases the electronic bandwidth and in turn 
enhances the relative strength of the Hubbard repulsion, 
bringing the material closer to a Mott insulating state. 
Recent experiments where K-(ET) 2 Cu[N(CN) 2 ]Br was re¬ 
versibly switched from a metallic to an insulating state 
by tuning the endgroup configurations are easily under¬ 
stood in our picture. In K-type materials that are not 
close to any phase transition, the effects of ethylene end- 
group disorder might not manifest as dramatically as in 
K-(ET) 2 Cu[N(CN) 2 ]Br. Based on our estimates of model 
parameters, k"-(ET) 2 Cu[N(CN) 2 ]C 1 might exhibit simi¬ 
lar behavior. 
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out to us and acknowledge support by the Deutsche 
Forschungsgemeinschaft through grant SFB/TR 49. Cal¬ 
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puting (CSC) in Frankfurt am Main, Germany. 
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